Intracavity weak nonlinear phase shifts with single photon driving. 
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We investigate a doubly resonant optical cavity containing a Kerr nonlinear medium that couples 
two modes by a cross phase modulation. One of these modes is driven by a single photon pulsed 
field, and the other mode is driven by a coherent state. We find an intrinsic phase noise mechanism 
for the cross phase shift on the coherent beam which can be attributed to the random emission times 
of single photons from the cavity. An application to a weak nonlinearity phase gate is discussed. 

PACS numbers: 



I. INTRODUCTION 



The field of quantum optics, now roughly four decades old, arose due to theoretical and technological innovations 
driven by the invention of the laser. From the very beginning single photon states were a central concern, especially in 
the context of photon antibunching. However it is probably fair to say that for most of those four decades the field has 
been largely concerned with Gaussian states, beginning with Glauber's work on coherence and subsequent work on 
squeezed states in the 1980s and 1990s. Much of the formalism of quantum optics, especially quasi-probability function 
methods, is based on systems driven by gaussian states. In the last decade however this has changed and there is now 
an emphasis on one and two photon sources Q, and their applications in quantum information processing 0, The 
study of single photon sources driving polarisable media is just beginning^ [f| 0, H, Q. 

We consider a Kerr nonlinear medium, inside an optical cavity, that couples two-modes by cross phase modulation, 
see figure [TJ One of these modes is driven by a pulsed field prepared in a photon number eigenstate, that is to say, 
each pulse has one and only one photon. The other mode is driven by a coherent state, which may be pulsed or a 
continuous wave. The objective is to compute the phase shift of the coherent field emitted from the cavity that is 
induced by the single photon. Schemes such as this form the basis of quantum gates based on weak nonlinearities fiol] . 
We also find an intrinsic phase noise, similar to that of [l~fl ] which can be attributed to the random emission times of 
single photons from the cavity. 
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FIG. 1: An optical cavity contains a Kerr nonlinear medium and is doubly resonant for two modes, labeled a,b. 
driven by a sequence of single photon pulses, while mode-b is driven by a coherent field. 



Mode-a is 



The cavity is doubly resonant for two modes, labeled a, b. The input to mode-a is a pulsed field at carrier frequency 
fl a with a single photon excitation per pulse. The input to mode-6 is a coherent field, either continuous or pulsed, 
at carrier frequency We first describe the single photon input states at the carrier frequency Q a . The cavity 
geometry ensures mode-a and mode-o are coupled to quasi-one dimensional modes outside the cavity. We assume 
all other modes are in the vacuum state and ignore them. In the continuum limit we define the positive frequency 
operator for the multi-mode field input to the cavity mode-a, at the input mirror, as 
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[ai(wi), a\(oj 2 )] = 5{loi - w 2 ) 
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Only some finite bandwidth, B of these modes are excited around the carrier frequency £! a . At optical frequencies, 
Cl a » B which enables us to set the lower limit of integration to minus infinity in Eq.|T]). In this form the moment 
n(t) = (aj(f)aj(t)) has units of s _1 . This moment determines the probability per unit time (the count rate) to count 
a photon at the input mirror between time t and t + di[l3j]. Likewise, we define the positive frequency input operators 
for mode-fe: 

h(t) = e -'^' / dubi{u>)e- iut (3) 

V 27T J-oo 



The multimode single photon state is defined by [13] 

/■oc 

|1) = / du v(w)at(w)\0) (4) 



Normalisation requires that dui\v(u;)\ 2 = 1 This state has zero average field amplitude but n(t) = |^(£)| 2 where 
u(t) is the Fourier transform of v(oj). For this state the function v(t) is periodic in the phase ut and it is not difficult 
to choose a form with a well defined pulse sequence. However care should be exercised in interpreting these pulses. 
They do not represent a sequence of pulses each with one photon, rather they represent a single photon coherently 
superposed over all pulses. Once a photon is counted in a particular pulse, the field is returned to the vacuum state. 
It would be more realistic to consider a field with one and only one photon per pulse, but for this paper the simpler 
version will suffice. A single photon state is a highly non-classical state with applications to quantum information 
processing. A review of current efforts to produce such a state may be found in 

II. QUANTUM STOCHASTIC DIFFERENTIAL EQUATIONS. 

We assume that the input field for mode-6 is a coherent state with a carrier frequency equal to fib and an amplitude 
of E(t). This is chosen so that the probability per time to count a single photon in this mode is |-E(i)| 2 for a unit 
efficiency detector. A simple way to model this is to include a classical driving term in the interaction Hamiltonian 
as Hd — e(t)(b + &'), with e(t) = s/KbEit) The input field to mode a is not a vacuum state, but non stationary pulsed 
single photon states. In this case an equivalent formulation in terms of the quantum stochastic differential equations 
enables us to refer directly to the fields external to the cavity. These take the formfl^, [3] 

^ = -i5 a a - ixhb{t)a - + y /K a ~a l (t) (5) 
-ie(t) - ixn a {t)b ~ ^-b + y^bi(t) (6) 



dt w A " v/ 2 

where the number operators are h a (t) = a^(t)a(t), fib = b^(t)b(t), and x is proportional to the third order optical 
nonlinearity, while S a = u a — Q a , where u> a is the cavity frequency for mode-a, and with K a , K n the photon number 
decay rates for modes a and b respectively. In this formulation the quantum noise terms in fact represent the input 
modes to the cavity. In this case, the noise term bi(t) is vacuum noise, however the noise for ai(t) requires careful 
treatment. We can also define output modes, in the interaction picture, 

a Q (t) = — = / duja (uj)e- wt (7) 



/2tt 

i r°° 

b Q (t) = —— \ dujb (uj)e- wt (8) 
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which represent the quantum noise terms for the time reversed versions of the quantum stochastic differential 
equations [3]. The input and output operators are related through the cavity modes 

a a (t) = y^a(t) - a,i(t) (9) 
b (t) = ^T b b(t) - h(t) (10) 



The standard approach to solving these equations is to take the Fourier transform and neglect the initial transients 
to obtain algebraic equations for a(w), a,i(u>) etc. If the equations are linear, the output statistics can then be directly 
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determined from the input statistics. However here the equations are non linear and this approach will not work. 
We will follow a different approach inspired by the quantum trajectory formulation in terms of realisations of the 
quantum stochastic differential equations. 

Our approach is based on replacing n a (t) by a classical stochastic process n a (t). This is a stochastic binary variable. 
The transition rate for the process n a = — * n a = 1 is the absorption rate for the single photon in the input state 
to mode-a. As discussed below, this is given by Eq.(fl7|). The transition rate for the process n a = 1 — * n a = is the 
emission rate for the single photon in the input state to mode-a. In effect this way of thinking treats the cavity for 
mode-a as if it was a photon detector. Of course no real record is available of this detection event, but as it is an 
irreversible process we can still think of it as a detection event in principle. 

We first calculate the absorption rate, that is, the probability per unit time that the cavity makes a transition 
from the vacuum state to a one photon state |0) a — > |l) a . The stochastic differential equation for mode-a is given 
by Eq(5). We will now make an approximation n b {t) rs n b (t) ~ |e(i)| 2 . This means we are ignoring the effect of 
quantum fluctuations of the coherent field, mode-b, on the single photon field. We will address the validity of this 
approximation in section HTT1 where we show that it is justified so long as \ is the smallest frequency in the problem. 
The stochastic differential equation for mode-a is then 

^ = -i X n b (t)a - ya + y^a^i) (11) 

where nj,(t) = K b \E(t)\ 2 is a known function. The formal solution is 

a(t) = a(0)e- K ^ 2 - 1 ^ + y^e-^*/ 2 -^*) f dt 1 'a^e^ 7 2 +«(*') (12) 

Jo 

with 

£(t) = X I dt'n b (t') + SJ (13) 



Given that the cavity field at time t — is not expected to be correlated with the input field at time t > 0, we see 
that 

d{n a ) 

—fa- = -Ka(n a (t)) (14) 

+ Ka e- K ° t / 2 +^ f d{{a\tf) ai {t))e Kat ' ,2 ~ im 
Jo 

+ Ka e- K " t / 2 -^) f ^'(a l t (t) ai (t'))e K " t '/ 2+l « t ') 
Jo 

The two time correlation function for the input field can be explicitly evaluated given the input state, Eq.f4|); 

(al(h)a i (t 2 ))=v*(t 1 )v(t 2 ) (15) 

The first term in Eq. (|14[) corresponds to emission of photons from the cavity at the Poisson rate K a . The second term 
is thus the absorption rate, or probability per unit time to absorb the single photon into the cavity, which we write 
as Pabs(t), so that 

PaUt) = K a e- Kat/2+m ! dt'v^t')^^'' 2 -^ + c.c (16) 



which is completely determined by the decay rate for mode-a, and the input fields for modes a and b. In the special 
case that x = and v(t) is real this reduces to 

Pabs (t) = 2 Ka e- K ^ 2 [ dt' V (t')u(t)e K ^'^ (17) 



If mode-b is driven by a continuous field, n b (t) can be replaced by its steady state value, n b . If we then detune the 
mode-a cavity so that 5 a — —X n b we can then set £(t) =0 and use the result in Eq.(fl7|) for the absorption rate in the 
presence of the nonlinearity. This means that the mode-a cavity is on resonance even when the Kerr nonlinearity due 
to mode-b is taken into account. 



4 



As an example we consider 



which corresponds to a single photon being emitted from a cavity, with damping constant 7, when prepared in a single 
photon state at time to, which we now set to zero, to = 0. The absorption rate is then given by 



PaUt) = (e-(--^ 2 - l) ( 7 e-^) (19) 

7 — K„ V / 



For both cases, n a > 7 and 7 > n a it is positive. 

If we substitute the absorption rate into the equation for the mean photon number in mode-a, Eq. (|14p . we can solve 
it directly to give 

(n a (t)) = , 47 "° e-"°' (l - e^-^) 2 (20) 

This is the ensemble average photon number in mode-a, averaged over all absorption and emission events. 

To assist with the interpretation of the various probabilities we first note that in our model we only calculate time 
dependent probabilities and rates; there is no spatial dependence. This is the point of view taken in the input-output 
formalism that has proved to be so effective in quantum optics. There are three mutually exclusive ways in which a 
photon can be detected between times t and t + dt: (i) it can be detected in the source with probability Pr s (t) , (ii) 
it can be detected in the cavity with probability Pr c (t) and (iii) it can be detected outside the source and the cavity 
with probability Pr a (t) = 1 — (Pr s (t) + Pr c (t)). In our model it is easy to see that 

Pr s (t) = e-T* (21) 
Pr c (t) = (h a (t)) (22) 

where the last equation follows from the fact that there can only ever be at most one photon in the cavity. It then 
follows that the rate of detection of photons outside the cavity must be 

Using Eq. (TT4")) together with the interpretation implicit in Eq (jXTJ) we then see that 

Po (t) = \v(t)\ 2 - Pabs (t) + n a (h a (t)) (24) 

On the other hand, the input-output formulation indicates that the rate of detection for photons after an interaction 
with the cavity is 

Po(t) = (n a (t)) = (<4(t)a (t)) (25) 

If we use the input /output relation in Eq.® we find that 

(n (t)) = \v{t)\ 2 - Pabs (t) + K a (n a (t)) (26) 

which agrees with Eq. (j2"4")l . This simply says that the rate of counting photons after the cavity is the rate at which 
photons are generated, minus the rate they are absorbed by the cavity plus the rate at which they are emitted from 
the cavity. In figure [5] we plot the source emission rate, the absorption rate and emission rate of the mode-a cavity 
for two values of n a . In both cases 7 = 1.0. 

We return now to the calculation of the phase shift on modc-b produced by the single mutual Kerr interaction with 
mode-a, driven by a single photon source. In a quantum trajectory picture, the number of photons in mode-a is a 
classical stochastic process n a (t) which takes the values zero or unity. The nonlinear quantum stochastic differential 
equation, Eq.©, can now be replaced by 

^=-ie{t)-ixn a {t)b-^-b + V^h(t) (27) 

We consider a continuous wave for mode-b and set the initial condition long before the single photon pulse reaches 
the cavity, so that mode-6 dynamics is given by Eq. (|27[) with n a = 0. In that case, the mode-6 will have reached a 
steady state which is in fact a coherent state with amplitude Poo — —j^- 
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FIG. 2: The plot of the source emission rate (solid), the absorption rate (dash-dot) and emission rate (dot) of the mode-a 
cavity for, (a) n a = 0.5, (b) n a — 0.9. In both cases 7 = 1.0 



The single photon input to modc-a will have entered the cavity at some random time and n a (t) = 1 for t>tA- 
At some Poisson distributed time, tg, the photon will be emitted from the cavity, at rate K a , so for t > t& we have 
n a (t) = 0. In that case, T — t^ — tj± = T, is a random variable corresponding to the dwell time of the mode-a photon 
in the cavity. Over this time interval the state of mode-6 evolves as a coherent state with amplitude, 

(3(T) = + , 4£X n e -(^/2+ix)T (28) 

where T is a classical random variable. We can determine the change in the phase of the field from the steady state 
value, Ap(T) = (3(T) - (3^ as 

A/?(T) = . 4e * f e -("»/2-Hx)r _ A (29) 

Hy ' K b (n b + 2i X ) V J V 1 

We expect that in typical situations that K b » x so that we can approximate 

A/3(T) = ( e —" T / 2 - 1) (30) 

This is the conditional change in the amplitude of the probe field from steady state, conditioned on a photon having 
entered mode-a. 

The statistics of T is determined by the absorption and emission of the photon by the cavity. Recalling that we 
are in a 'quantum trajectory' picture in which n a (t) is zero or unity, the conditional emission probability given that 
n a (t) = 1 is a Poisson process with rate n a - The probability distribution for T is then given by 

p(T) = n a e-^ T (31) 

The average dwell time, T, is then simply (K a ) _1 . The change in the amplitude, Eq. (f3"0"l) will thus be greatest when 
Ka/Kb << 1< Furthermore the fluctuations in amplitude change will also be small. We can test this intuition by 
simulating the Poisson process for T. In figure [3] we show histograms of \A/3(T)\ given in Eq. (|3"D]) . with t\j ti\ = 1. 
In figure[3Ja) we show an example in which K b / n a < 1. The wide variance in the amplitude shift is evident, whereas 

in figure Efb), for which K b /n a > 1 and x/ K b, the change in amplitude is peaked at -4£. The fluctuations in the 

K b 

displacement A/3(T) result from the photon leaving mode-a before mode-b has a significant chance to respond. 

We can now compute the ensemble average of the conditional amplitude change, A/3(T), over the distribution of 
waiting times, T . The mean is given by 

{Kb T ) 

where B = In the limit n b >> n a we find that the conditional mean is —B and the variance goes to zero. This 

K >> 

result is a conditional mean conditioned on the fact that the single photon does indeed enter the cavity. We now 
wish to compute the unconditional average of the change in the amplitude for the probe, at time t, averaged over all 
emission and absorption events for mode-a. We will do this using the cascaded systems formalism of quantum optics. 



while the variance is 
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FIG. 3: A histogram of the changes in the amplitude of mode-b from the steady state value for many realisations of the time, 
T, over which the single photon resides in the cavity for, 4ex/nl = 1 and (a)K(>/Ka = 4.0, (b)Kb/K a = 0.25. The number of 
trials is 50, 000. 



III. CASCADED SYSTEMS APPROACH. 



We can also describe the system in terms of an ensemble average over the random absor ptio n and emission times of 
the single photon using the cascaded systems method of Carmichael and also Gardiner [l5lll6|| . One cavity, labeled by 
annihilation operator, c, is called the source, the output from which is coupled uni-directionally into the second cavity 
for the two modes a and b coupled by a cross Kerr interaction as before. The system is described, in the interaction 
picture, by the master equation 



dp 
~di 



—iS a [a'a, p] — ie[b + b' , p] — ixia^ab^b, p] 

+jV[c]p + K a V[a]p + K b T>[b]p + y/~/K a ([cp, a 1 ] + [a, pc 1 ]) 



(34) 



We assume that the source cavity starts at time t = in a one photon state |l) c and first consider the situation with 
X = 0, i.e we only consider mode-a. The equation of motion for the number in the driven cavity is then 



d(n a ) 
dt 



= -K a (n a ) - V7'« a ((a t c) + (c f a)) 



(35) 



can be solved directly and substituted into Eq. (|35l) to obtain 



The second term is the absorption probability in analogy to the second term in Eq (|14p . The equation of motion for 



(36) 



d(n a ) 
dt 



(7 - «o) 



( 7 +«„)t/2 / : _ e _( 7 -« a ) t /2) 



This implies the absorption probability per unit time for a single photon to be absorbed from the source cavity is the 
second term and this agrees with the result in Eq. (fl9| . We can solve Eq. (f36|) directly to obtain the result quoted in 
Eq.®. 

We assume that mode-b starts at t = in the steady state coherent state |/3oo) while mode-c starts in a single 
photon state |l) c while mode-a starts in the vacuum state. We now make a canonical transformation to a displacement 
picture p{t) — D(f3 OQ )p(t)D' i (Poo)- This shifts mode-b to a new equilibrium position so that at t = mode-b starts in 
the vacuum state as well as mode-a. The master equation in the displaced picture becomes, 



= -ixW a-iPootf + (3^b, p] - ixi^arfb, p] 
+~fV[c]p + K a V[a]p + K b V[b]p - 



(37) 



[a,pc f ]) 



We also assume that the detuning 5 a has been chosen to cancel any systematic detuning due to the back action of the 
probe field. We remind the reader at this point that using the master equation in Eq. ([38|) describes the unconditional 
dynamics of the system, averaged over all possible emission and absorption events from each cavity. We thus do not 
expect to be able to reproduce the conditional phase shifts discussed in the previous section. However we can calculate 
ensemble average phase shifts on the mode-b using this master equation. 
The equations of motion for the mean amplitude in mode-b is 



d(b) d 
dt 



(38) 
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where the d subscript reminds us that we are in a displaced frame in which [3^ has been subtracted. The equation of 
motion for (n a ) is still given by Eq.p5p. however there is now an implicit \ dependence, as the equation of motion 
for (c'a) depends on \- If we neglect this dependence and approximate (n a ) in Eq. ((38|) with the solution in Eq. ([20|) 
we are in effect seeking a solution to linear order in \. We also assume that \ftoa\ >> 1 and neglect the final term 
in Eq. (|38p . This is justified as the state of mode-a is at most one photon away from the vacuum while to lowest 
order in \, mode-b is near the vacuum state in the displacement picture. Neglecting this term means we are ignoring 
correlation between mode-a and mode-b and is consistent with the assumption that mode-b remans in a coherent state 
as we assumed in the section [TXJ The solution for (b) in the displacement picture is then easily obtained. However in 
section [H] we saw that best performance is likely to require that mode-b is slaved to mode-a, that is to say Kb » n a . 
In that case we see that the change in amplitude of mode-b in the ensemble average is, to linear order in x, 

Ap(t) = (b)d = -^f(n a (t)) (39) 
K b 

Of course this eventually decays to zero as the photon in mode-a eventually leaves the cavity and mode-b simply 
relaxes back to the steady state coherent state of amplitude /?oo- In contrast the condition mean in Eq. ([32|l assumes 
that the photon has entered the cavity but has not left the cavity. 



IV. DISCUSSION AND CONCLUSION. 



One of the important applications of the model of this paper is to enable a two-qubit parity gate based on an 
intracavity weak nonlinearity, (see figure 2]). The qubits are encoded as dual rail single photon states, for example, 
each rail could correspond to orthogonal polarisation states of a single spatial mode. In figure [4] we distinguish the 
modes of each rail using an overbar notation, ak,dk, where k = 1,2 for each qubit. Only one rail of each qubit is 
directed towards the cavity inside of which we model the interaction as 

H = hx{a\ax - 4a 2 )6 t 6 (40) 

We are not suggesting this as a realistic description of a mutual Kerr nonlinearity but rather as a minimal model for 
a two qubit parity gate for dual rail single photonic qubits. Our objective here is to assess the effect on the gate of 
the random interactions times of each photon with the coherent field in the cavity. The parity gate works as follows. 
If the number of photons in each mode are the same, that is to say both or both 1, there can be no phase shift 
on the coherent probe field. If however the photon number in each mode is different, there is an equal and opposite 
phase shift for the two cases n% = l,ri2 = Q and n% = 0,n2 = 1. This assumes that each photon interacts with the 
probe field for the same time. However as we saw in the single mode case, the absorption and emission of photons 
from the cavity are stochastic events. Furthermore in this model the stochastic events for each mode are independent 
(if we ignore the back action of the probe field on the single photon modes). 




FIG. 4: A model of a two-qubit parity gate for dual rail code with single photon states. The qubits are labeled as Q^. The 
probe field is a strong coherent state. The box represents a cavity resonant with both the probe and the qubit modes. One rail 
of each qubit does not pass though the cavity. 

Using the results in the previous section, the probe field is injected as a coherent state with amplitude t{t)/ ^/kZ, in 
which case it will remain a coherent state with amplitude (3{t) given by 

^ = -ie(t) - » X (m(t) - n 2 {t))P - ^9 (41) 



8 



where n%(t), n,2{t) are well approximated by independent stochastic variables that take values 0, 1 and correspond to 
the number of photons in cavity modes a\,a2 respectively. We again assume that e(t) — e is a constant and that the 
probe field has relaxed to a steady state long before the qubit photons enter the cavity. We also assume that the 
qubit photons are identical synchronous single photon pulsed states given by Eq. 2J As before we expect the best 
performance to occur when the probe field is heavily damped so that it is slaved to the qubit modes. In the typical 
case for which Kb > > % we then find that the change in the amplitude of mode-b from the steady state value is given 
by 

A/3(t) = -lf K(t)-n 2 (t)) (42) 

For a perfect parity gate we would like this amplitude change to be zero for the two cases in which both modes had 
either zero or one photon. However it is now clear that this can only be true on average. From shot to shot, we are 
not guaranteed that the stochastic variables both take the same values even when both qubits have zero or one photon 
in both of the cavity rails and so, from shot to shot, the amplitude change will fluctuate around the ideal case of 
zero. Likewise for the case in which the number of photons in each of the cavity rails is different. In that case would 
like there to be a significant equal and opposite amplitude change, but now there will be fluctuations superimposed. 
A detailed error analysis of quantum circuits based on the effect of these fluctuations will be the subject of a future 
paper. 

In this paper we have presented two approaches to nonlinear intracavity systems driven by single photon states. 
One approach is based on a generalization of the quantum Langevin approach to the case where the input noise 
describes single photon excitations of the extra-cavity field. The second approach is based on the cascaded system 
master equation of Carmichael[15J and Gardiner[16]. In both cases we see intrinsic noise arising from the random 
absorption and emission times of the photon. These processes are however modeled as coherent interactions between 
the fields inside and outside the cavity. As a simple application of the formalism we discussed a two-photon parity 
gate based on weak nonlincarities, however many more examples could be given. Of particular interest is the case of 
single photons interacting with cavities containing single atomic systems. This situation arises very naturally in the 
context of quantum repeaters [13], a key component of a quantum internet based on single photonics [l8j]. Given the 
effort going into making controllable single photon sources we expect that single photon quantum optical systems of 
the kind discussed in this paper will play an increasingly important part in the future of the field. 
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